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Abstract 

The very small scales of isotropic, Navier-Stokes turbulence at Reynolds 
number 1Z\ ~ 15 are studied by high-resolution direct numerical simulation 
(DNS) and by integration of the direct-interaction (DIA) equations. The 
DNS follows the tail of the energy spectrum over more than thirty decades 
of magnitude. The energy spectrum in the far-dissipation range 5kd < k < 
Wkd is well-fitted by k a exp(— ck/kd), where kd is the Kolmogorov dissipation 
wavenumber, a ~ 3.3 and c ~ 7.1. For values of m that emphasize the far- 
dissipation range, the fields (— V 2 ) m u exhibit strong spatial intermittency, 
associated with gentle spatial variations of the lower-/c part of the velocity 
field. DIA analysis gives a prefactor A; 3 and an exponential decay more rapid 
than DNS. Averaging over an ensemble of DIA solutions, suggested by the 
observed intermittency, removes some of the discrepancy. 
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The smallest scales of incompressible, isotropic Navier-Stokes turbulence are associated 
with the far-dissipation range of wavenumbers k 3> k d , where kd = e 1 / 4 ^ -3 / 4 is the Kol- 
mogorov dissipation wavenumber, e is the rate of dissipation of hydrodynamic kinetic energy 
per unit mass and v is kinematic viscosity. Most of the dissipation takes place at k < kd- 
The wavenumbers k 3> kd have attracted attention for a number of years. There has been 
controversy concerning the asymptotic form of the energy spectrum as k — > oo. The small- 
est scales are of further interest because they display strong intermittency even at Reynolds 
numbers so low that there is no basis for a fractal cascade. 

It is reasonable to assume that the wavenumbers k^> kd represent spectral tails of flow 
structures of spatial scale > 1/kd- An analogy is the exponential spectral tail of a shock that 
obeys Burgers' equation. There are other possibilities. One is that very-high-wavenumber 
excitation comes mostly from exceptionally strongly strained regions that give rise to ob- 
served exponential-like skirts of the probability distribution function (pdf) of vorticity. 

A number of authors have discussed kinetic energy spectra for k 3> kd of the form 

E(k) oc f(k/k d )exp[-c(k/k d ) n ] (1) 

where c is a constant, / is a weak function of k/k d and 1 < n < 2 The direct 

interaction approximation (DIA), a perturbative treatment, gives n = 1 and / oc (k/kd) 3 
@. Perturbation approximation can be justified for k kd because the mode amplitudes are 
very nearly the linear response, under molecular viscosity, to quadratic forcing by modes of 
lower k ||. The additional, unjustified assumption in the DIA analysis is that the statistics 
of wave-vector amplitudes for k ^> kd are nearly Gaussian. The plausible effect of the 
intermittency actually present in the far-dissipation range is to enhance mean nonlinear 
transfer and thereby raise the spectrum level above that predicted by DIA. 

Foias, Manley and Sirovich || have shown that n > 1, under certain assumptions of 
smoothness of the velocity field in a finite box. In view of this inequality, The DIA results 
suggest that n — 1 is exact for a finite box. The particular form 

E(k) oc k a exp(-ck/k d ) (2) 

seems consistent with a body of experimental and computer data 0. 

Strong intermittency in the far-dissipation range at modest Reynolds numbers was pre- 
dicted some years ago on the basis of a simple physical argument [§]: E(k) falls off steeply 
for k ^> kd- Consequently, fluctuation, on spatial macroscales, of parameters like kd in ([!]) 
yields spatial intermittency at scales 0(1/ k) that increases without limit as k/kd — ► oo. 

An ongoing computation project has achieved resolutions up to 512 3 (wavenumber range 
1 < k < 256) on a CM-200 computer at Los Alamos National Laboratory ||. The direct 
numerical simulation (DNS) described in this paper is limited by arithmetic precision rather 
than its resolution of 256 3 (wavenumber range 1 < k < 128). A nominal steady state was 
maintained by forcing confined to k < 3, at a level determined to give the desired value 
of the Taylor microscale Reynolds number TZ\. The Taylor microscale A for an isotropic 
turbulent flow is a length defined by A = (lSz/Vg/e) 1 / 2 , where t> 2 is the mean-square velocity 
in any direction, and 7Z\ = v^X/v. 

The solid line in Fig. [1] represents the time-averaged wavenumber spectrum of kinetic 
energy E(k) for a run with v = 0.026, 7Z\ ~ 14.9, k d ~ 9.65. The very close fit of the 
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high- A; part of the spectrum to a straight line (exponential decay) is apparent. Fig. |] shows 
kdlnE(k)/dk vs k for this data. If E(k) has the form (g), this plot is a straight line whose 
slope is —c/kd and whose intercept on the vertical axis is a. The straight line in Fig. |2] is a 
least-squares fit to the data over the range 50 < k < 100; it gives c ~ 7.1, a ~ 3.3. Equation 
(0) seems well supported. The confidence level for a is not high, because it is not certain 
that the wavenumber range is long enough to give strictly asymptotic results. It cannot 
even be asserted that f(k/kd) is exactly a power. However, the value n = 1 in (ffl) does seem 
strongly favored. 

A previous study J7J gave a negative value for a. The range used for fitting in J7J was 
0.5kd < k < 3kd, which is too low to give asymptotic behavior. Note that in Fig. |2| the 
data points curve downward at small k. An intercept a = corresponds to tangency at 



k = 30 « 3k d p 



The dotted lines in Fig. [I] show the spectra in three subregions of the cyclic box, defined 
by the x-space filter exp[— |x — x c | 2 /(L/32) 2 ], where L = 2tt is the box size and x c is the 
subregion center. The nominal linear dimension of a subregion is thus L/16. In order to 
sufficiently reduce errors from chopping, the cyclic box was repeated in each direction to 
give a total of 3 3 replicas before the filtering. The filtered field was transformed to fc-space 
and subjected to solenoidal projection before the spectrum was computed. An effect of the 
filtering is marked depression of spectrum level for k ^ 16. 

We have found that the distribution of spectral slopes over a set of 64 subregions, evenly 
spaced in the cyclic box, is consistent with the picture of intermittency || in which the 
parameters describing the spectral tail are slowly-varying functions of spatial position. The 
three subregion spectra plotted in Fig. [I] are those with minimum, median and maximum 
values of E(k = 60). Note that the maximum and minimum values of E(k = 60) differ by 
a factor of over 10 10 , despite the relatively small difference in slope. 

Fig. H] shows pdf's of the differentiated velocity fields (— V 2 ) m Mj, averaged over the three 
components i = 1,2,3. Note the marked increase of intermittency with m. Fig. ^ is a 
visualization of the field V 8 «i. The shaded surfaces are where the absolute value of field 
amplitude equals twice its root-mean-square value. The evident three-dimensionality of these 
regions suggests intermittency that is associated with gentle spatial variation on principal 
dissipation scales rather than with exceptional regions that are strongly strained into thin 
sheets or tubes. The flatness ((V 8 wi) 4 )/ ((V 8 wi) 2 ) 2 is 57. This large value is associated with 
the sharp peak of the pdf at zero amplitude, rather than with the broad skirts at large 
amplitude values. The latter represent probabilities too small to affect low-order statistics. 

The spectral support of the field V 8 «i is effectively confined to the range 15 < k < 40. 
A field with spectral support confined to 50 < k < 100, the region where E(k) is accurately 
proportional to k a exp(—ck/kd), can be constructed by applying the filter exp[— {k— 75) 2 /200] 
to the wave- vector transform of wi(x) and transforming back to x space. Fig. [5] shows the 
regions where the absolute value of the amplitude of this field exceeds twice its root-mean- 
square value. The field V 24 «i has approximately the same spectral support, but it is too 
noisy to give clean visualizations. 

The similarity between Figs. (| and [| is striking. The regions of high intensity are in the 
same locations and have similar shape, but are smaller in Fig. |5]. This behavior suggests that 
the spectral support of the very small scales represents the spectral tail of larger structures. 
Similar behavior is exhibited by repeated differentiation of the velocity profile of a Burgers 
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shock. 

We have integrated the DIA equations with the same forcing and viscosity as in the DNS. 
The dashed line in Fig. [I] shows the result for E(k). The DIA spectrum tail falls within the 
range of DNS subregion values, but below the median. There are two obvious causes for 
discrepancy between DIA and DNS. One is the depression of high-fc energy transfer in DIA by 
sweeping effects in response functions Q. At k — 75, the sweeping decorrelation frequency 
Vok is about 30% of the viscous decay frequency vk 2 . Two-point closures that are invariant 
to random Galilean transformation do not display the strong sweeping decorrelation at large 

k H- 

The second cause is the observed intermittency, which is not captured by DIA. Consider 
the mean spectrum obtained by averaging over a finite-size ensemble of DIA solutions, 
or other closure solutions, with different values of spectrum level in the forcing region. As 
k — > oo, this mean spectrum is dominated by the ensemble member with the smallest value of 
c/kd- Thus the effective value of c/kj is decreased by the averaging. The prefactor exponent 
a = 3 is unchanged by the averaging. This latter value is a property of a broad family of 
two-point closures 0. The averaging over closure predictions is not explored further here, 
but it remains possible that simple statistical distributions of the input parameters to the 
closure can give satisfactory approximations both to E(k) and to the large flatness factors 
of the fields (-V 2 ) m u. 

The procedure of averaging over solutions with different parameter values has an ana- 
lytical implication that does not depend on closure approximation. Suppose that the box 
containing the flow is of infinite size, but with spectral support confined to wavenumbers 
above some wavenumber ko. A local wave- vector analysis centered on an arbitrary point 
x c can be carried out by applying the filter exp(— /c 2 |x — x c | 2 ) to the velocity field and 
performing solenoidal projection on the result. 

Suppose that the spectral tail thereby associated with a region of the box with dimensions 
0(1/ ko) has the asymptotic behavior (g), where k d is now a weak function of x c . If any value 
of fcrf(xc), however large, occurs for some x c , then averaging over the entire infinite box will 
give something slower than exponential decay of the spectrum as k — > oo. For example, if 
the pdf of k d over the box were oc exp(— k d /k 2 ) as k d — > oo, where k c is a constant parameter, 
then averaging of (fj) over this pdf would give a full-box spectrum E{k) oc exp[— (2ck/k c ) 2 ^ 3 ] 
as k — > oo, apart from an algebraic prefactor. Finally, if the size of the box were made finite 
but 3> l/ko, then this full-box spectrum would change over from stretched-exponential to 
simple-exponential decay at some wavenumber 3> k c . The bound obtained by Foias, Manley 
and Sirovich || assumes a finite box. 
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valuable discussions with S. Kida, O. Manley, R. S. Rogallo, F. Waleffe, and Y. Zhou. D. W. 
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the Department of Energy, the National Science Foundation, the Defense Advanced Research 
Projects Agency, and the Office of Naval Research. The computations were performed at 
the Advanced Computing Laboratory, Los Alamos National Laboratory and at the National 
Center for Atmospheric Research. 
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FIGURE CAPTIONS 

FIG. 1. Linear-log plot of E{k) vs. k. Solid line, DNS spectrum; dotted lines, subregional 
spectra with largest, median and smallest values of E(k = 60); dashed line, DIA spectrum. 

FIG. 2. Lhe function kdlnE(k)/dk vs k for the DNS spectrum. The straight line is a 
least-squares fit to the data points for 50 < k < 100. 

FIG. 3. Pdf P(w) of the field w = (-V 2 ) m u; for m = (dashed), m = 2 (dotted) and m = 4 
(solid), averaged over i = 1, 2, 3. 

FIG. 4. Perspective view of the surface where |V 8 «i| equals twice its root-mean-square value. 

FIG. 5. Perspective view of the surface where the absolute value of the filtered field with 
spectral support centered at k = 75 equals twice its root-mean-square value. 
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